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The opening of the meson factories twenty years ago provided nuclear physics with new 
beams, higher momentum transfers, and new opportunities for precision measurements. The 
resulting changes in nuclear physics were substantial, altering not only the range of physics 
issues identified with the field but also the manner and size of the collaborations that do nu- 
clear physics. Inspired by the talks of this symposium, I discuss some of the accomplishments 
as well as some of the goals not yet reached. 



1 Introduction 

It is a great pleasure to take part in this symposium, a celebration of twenty years of 
new physics from LAMPF and the other meson physics facilities. As in the cases of others 
who have spoken, LAMPF is a special place for me. My first workshop talk was given here, 
shortly after I joined T5 as a postdoc early in 1977. It was an occasion I remember fondly, as 
the few slides I showed stimulated a lively debate between Gerry Brown and Arthur Kerman. 
The audience seemed to give me some of the credit for helping to provoke their entertaining 
exchanges. Earlier my thesis, an investigation of inclusive nuclear interactions of neutrinos 
from stopped pions, had been inspired by Ray Davis's proposal to calibrate his solar neutrino 
detector at the LAMPF beamstop. 

Despite the assigned title for my talk, I hope the audience is not expecting a summary: 
the excellent talks of the past two days require none. I would instead like to make a few 
remarks that were inspired by what I have heard at this symposium, but which also connect 
to my own research. So, in a sense, this is more of a personal tribute to the meson facilities 
for the intellectual stimulation they have provided and continue to provide. 



2 Pion-nucleus physics and nuclear structure 

The experimental achievements of the meson factory 7r-nucleon and vr-nucleus programs 
were summarized in the excellent talks of Bill Gibbs, Jerry Miller, Bob Redwine, and Anna 
Hayes. They described a number of substantial achievements: the general accord achieved on 
low and medium energy vr-p data; credible results on isospin tests; discovery of the relative 
transparency of the nucleus to pions with E ~ 50 MeV; the successful program to exploit the 
cr(7r"'")/o"(7r^) selectivity to probe the isospin character of nuclear transitions; and mapping 
the isovector monopole resonance. 

The Bethe Panel, whose endorsement of the physics potential of meson facilities was 
crucial to the creation of LAMPF, was quite optimistic about the promise of pions as a 
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nuclear structure probe, especially in their sensitivity to less well understood aspects such 
as neutron distributions. It is clear that much progress toward this goal was made in the 
past twenty years: the elementary cross sections were measured quite well; theory produced 
optical potentials that accurately described elastic scattering; and a substantial body of 
data was gathered on single and double charge exchange reactions leading to specific nuclear 
states. Yet, despite this progress, the pion appears not to have furthered our understanding 
of nuclear structure to the extent envisioned by the Bethe Panel. 

While it is easy to attribute this situation to the more complicated reaction mechanism 
of strongly interacting probes, the parallel nuclear structure studies at electron facilities like 
Bates also produced puzzling systematics. For example. Ml form factors measured above q 
~ 400 MeV [1] tend to remain well above theoretical predictions, even though the electron 
probe is presumably perturbative and the nuclei in question, such as those in the Ip shell, 
can be treated with microscopic methods like the shell model. 

A common aspect of these Bates measurements and much of the pion physics program 
at LAMPF is that we are probing the nucleus at momentum transfers large compared to the 
Fermi momentum. I would like to argue here that some of our difficulties in understanding 
the nuclear response at high q^ could be attributed to associated shortcuts we make in theory. 

While the shell model is often viewed as a phenomenological technique, it is important 
to remember that it is at least qualitatively motivated by more rigorous many-body theory. 
Quite generally a problem in a large (or infinite) Hilbert space can be equivalently described 
in a smaller model space, provided that the interaction used in the model space takes into 
account the degrees of freedom that were "integrated out [2]." The correct Vcg-cctivc would 
reproduce the exact eigenvalues and eigenfunctions (or rather their restrictions to the model 
space) in the smaller model-space calculation. Now the full calculation of Veffective is a task 
fully 

equivalent in complexity to the original A-body problem. The hope in nuclear physics is 
that an expansion of Vcgcctive in clusters, or equivalently powers of the nuclear density, might 
rapidly converge [3] . If the convergence is sufficiently rapid, one might then have a tractable 
strategy for solving the A-body problem. 

Unfortunately the shell model has at best a tenuous connection to such rigorous many- 
body efforts. Most often the shell model potential is simply fitted to spectroscopic data. 
But lacking a diagrammatic basis for the effective interaction, one has no procedure for 
making the corresponding corrections that generate the appropriate effective operators from 
bare operators, a necessary step in evaluating any observable. In effect, the model has lost 
most of its predictive power. The alternative approach, calculating the effective interaction 
directly from the NN interaction, is sometimes carried out, but generally only in lowest 
order, summing the two-nucleon ladder diagram, or G-matrix. Often phenomenological 
adjustments are made to the results for the purpose of improving agreement with experiment. 
Matrix element evaluations with consistent treatments of effective interactions and operators 
are almost never attempted. 

Part of the preference for phenomenology may have its roots in the early years of the 
shell model, when numerical techniques were severely limited by the available computers. 
For example, convergence problems associated with intruder states [4] were consider a major 
obstacle; yet modern workstations can easily handle large-basis shell model calculations with 
many degrees of freedom, such as a full Ahui calculation of ^^O [5], thereby accounting for 
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all of the principle excitations affecting the low-lying spectrum. Furthermore, powerful new 
techniques could be applied to the problem of calculating higher order contributions to the 

effective interaction. 

Motivated in part by an INT program organized by Bruce Barrett and James Vary, two of 
us in Seattle are trying our hand in this old game. The starting point is the Bloch-Horowitz 
equation [6] for the effective Hamiltonian 

/^effective = (T + + (T + T/)-^_^l^g(T + V) (l) 

where T is the relative kinetic energy operator, V is the NN potential, and Q projects 
onto the excluded, high-momentum space. The hope in nuclear physics, which dates back 
to arguments made by Brueckner [3], is that the second term in this equation might be 
expandable in a series in the nuclear density, corresponding to the ladder diagrams for two, 
three, four, nucleons interacting at short range. 

The Bloch-Horowitz equation looks very difficult to solve because the energy E appearing 
in the Green's function is the exact eigenvalue. Thus, for each desired state, one would have 
to solve this equation selfconsistently, presumably iterating until convergence is achieved. 
If one had to reevaluate the effective interaction for each trial E, the algorithm would be 
clearly quite inefficient. 

However the Lanczos algorithm, quite commonly used to find the extremum eigenvalues 
and eigenvectors of very large matrices, has a less well known application to linear operator 
inversion [7] that circumvents this difficulty. The Lanczos algorithm solves the moments 
problem for the distribution of H over its eigenvalues, determining the 2n — l moments of H 
after n iterations. These moments can be determined from the tridiagonal Lanczos matrix. 
From knowledge of these moments one can calculate the effects of operators like the Green's 
function, which sums over the full spectrum. One finds 

^ ■vi)c^gi{E)\v,) + ... + g^{E)\vr,) (2) 



E-Q{T + V) 

where the \vi) are the Lanczos vectors and the gi are continued fractions determined by 
E and by the entries in the Lanczos matrix. If one chooses the starting Lanczos vector 
to be Q{T + V)\I) where |/) is a state in the model (included) space, this expansion 
can be used to determine the matrix elements of i^eflfective- Because the energy E is merely 
a parameter in the continued fractions, once one carries out the Lanczos expansion, the 
effective interaction is known as a function of E. No additional computer time is required 
to iterate on E. As the Lanczos expansion itself is quite efficient in determining the Green's 
function, typically converging in about 50 iterations, this becomes a very powerful method 
for solving the Bloch-Horowitz equation. 

We believe it may be possible to calculate the effective interaction in a cluster expansion 
through the three- and, possibly, four-nucleon ladder, thereby correcting the usually two- 
nucleon G- matrix through second order in the nuclear density. It is an open question whether 
this is sufficient to achieve convergence. 

This sketch of how one might attempt to generate i^eflfective given modern workstations 
and efficient algorithms, could be repeated for effective operators. This is the reason for 
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including the discussion here. Returning to the example of exclusive nuclear form factors 
measured at Bates, one can envision operating on a model-space wave function by the mag- 
netic operator Ml(g^) transferring, for example, twice the Fermi momentum. The resulting 
amplitude has a very small overlap with the model space relative-coordinates wave function 
describing the final state, as this wave function contains only long wavelength components. 
The scattered nucleon must interact with one or more neighboring nucleons, sharing the 
momentum, before the overlap becomes appreciable. That is, it is the effective contributions 
to the operator that should dominate the response. Therefore, if we fail to correct both 
our interactions and operators consistently, we will lack the basic ingredients for a realistic 
calculation of nuclear responses at high g^. And it is clear, as increases, we will need to 
include increasingly complicated correlations in our effective operator to achieve a realistic 
result. The arguments here are very much like those that lead to quark counting rules for 
form factors at asymptotically large q^. 

I believe our lack of effort on this problem - calculating consistent effective operators and 
interactions from realistic NN interactions - may have limited the meson factory (and other 
medium energy facility) impact on our understanding of nuclear structure. 

3 Neutron Physics 

Stuart Preedman, David Bowman, John Browne, and Steve Vigdor described many of 
the elegant symmetry tests that have been done with low and medium energy neutrons, 
as well as the prospects for future experiments at LANSCE and LISS. The results include 
remarkably precise limits on CP- and T-violation from neutron electric dipole moment and 
D-coefficient /5 decay measurements, as well as the demonstration that parity-nonconserving 
resonance mixing can probe chaotic properties of the nucleus. 

One interesting issue not discussed is the relation of neutron and atomic electric dipole 
moments measurements. The neutron edm limit from Genoble [8] and Gatchina [9] 



As the ^^^Hg atom is probed with electric fields of about 10'* volts/cm, this corresponds to a 
change in the ground state energy of the atom of about 10~^^eV, which is approximately the 
interaction energy of two electrons separated by one light-month! This certainly illustrates 
the extraordinary precision of the limits above. 

Although the ^^^Hg edm limit is about sixty times smaller than that of the neutron, 
deciding which is the more stringent test of underlying sources of CP violation is quite 
complicated. The nuclear edm of an odd-A nucleus is generated both by the edm of the 
unpaired valence nucleon and by a nuclear polarizability: a CP-odd NN interaction produces 
a parity admixture in the ground state that allows absorption of a CI (J=l multipole of the 
Coulomb interaction) photon. The latter is generally the dominant contribution to the 
nuclear edm [11]. Depending on the nature of the underlying CP- violating Lagrangian, the 



^8 • 10~^*'e cm 



(3) 



can be compared to the Seattle group's [10] limit on the mercury edm 




(4) 
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entire nucleus may contribute to the polarizing nuclear mean field, or the effects may only 
arise from the unpaired neutrons. This means that the question of relative sensitivity of 
neutron vs. atomic limits is quite dependent on the CP violation model. 

Further complicating the comparison is the Schiff shielding [12] of the atom: for a point 
nucleus in a neutral atom, the atomic interaction energy linear in the edm vanishes. The 
atomic cloud sags in the applied field to cancel the net field at the nucleus. (If this did not 
happen, the nucleus of a neutral atom would be accelerated by the external field.) In a real 
atom a nonzero interaction is found because of the nuclear finite size, hyperfine effects, and 
relativistic corrections. The residual effects grow rapidly with Z, making heavy atoms the 
favorites of experimentalists. 

The net result - nuclear coherence in the polarizability, the loss of sensitivity due to 
electronic shielding, and the factor of 60 smaller ^^^Hg edm limit - is comparable sensitivity 
to CP-violation sources that are coherent over the entire nucleus (such as the KM phase), 
and a modest advantage in favor of the neutron edm when the coherence is only over valence 
neutrons (such as in the case of the QCD 9 parameter) [13]. 

This situation may soon change. Recent improvements in the nuclear measurements 
have been very rapid, outpacing the corresponding improvements in the neutron edm limit. 
Furthermore, there are several interesting nuclei that should have greatly enhanced CI po- 
larizabilities, in one case by an estimated factor of 10^ [11]. While the best cases are not 
attractive experimental candidates, it is quite possible that this could change as new tech- 
niques (e.g., trapping experiments) [14] are developed. It is exciting that such advances 
might allow experimentalists to probe edms near the 10"^-^ e cm range predicted to arise 
from the standard model KM phase. 

4 Neutrinos 

Peter Nemethy, Peter Rosen, and Michael Shaevitz described the marvelous LAMPF 
contributions to neutrino physics (E31, E225, E645, and LSND). This continues an even 
longer tradition at Los Alamos that includes the discovery of the neutrino by Cowan and 
Reines, for which Reines won the Nobel Prize this year. The LAMPF experimental program 
provided crucial tests of additive vs. multiplicative conservation laws for lepton number and 
constrained the sign of the charge current - neutral current interference in (z/, e) scatter- 
ing. Most exciting is the current experiment, LSND, with 52 events attributed to neutrino 
oscillations [15]. 

One reason for the intense interest in neutrino physics has to do with the puzzling pattern 
of masses in the standard model. If this model is part of some grander, more unified scheme, 
one would anticipate finding larger multiplets in that extended scheme, perhaps uniting 
particles by generation, e.g., (e, i/g, <i). But as these particles presumably would have 
similar couplings to the fields generating masses, the much hghter mass of the Ue then seems 
a puzzle. However us are distinguished from the other particles in the multiplet because 
they have no additive quantum numbers that must change sign under particle-antiparticle 
conjugation. Thus in addition to the Dirac masses of the other fermions, neutrinos can have 
Majorana masses. In the popular seesaw mechanism, the introduction of a large right-handed 
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Major ana mass Mr leads to light and heavy neutrinos 



niD 



(5) 



rriL ~ TTiD 



niH ~ Mr 



(6) 



with the familiar z/g primarily identified with the light eigenstate rriL. Thus one sees that 
small neutrino masses are explained by a parameter, ttid/Mr, reflecting familiar and new 
mass scales. It follows that the ability to probe small neutrino masses becomes a search 
for very high scales Mr characterizing new physics. If we take Mr ~ 10^^ GeV, a typical 
grand unification mass, and ~ m^p quark, one finds a third generation neutrino mass of 
rriu^ ~ 10~^ eV, a number accessible to neutrino oscillation experiments. 

I think it would be a great disappointment if the LAMPF neutrino program were to come 
to an premature end: the beamstop remains our most intense terrestrial source of low-energy 
UgS, i/^s, and P^s. In the near term, the moving of LSND ought to be given very high priority: 
the associated cost of $5M is quite modest on the scale of the average neutrino experiment. 
The current results of that experiment make it imperative that every step be taken to test the 
results. On the longer term, I feel it is very important to the nation's nuclear physics program 
to keep a neutrino presence on the mesa. The neutron physics material science program at 
LANSCE calls for several possible upgrades. There is no reason that neutron and neutrino 
programs should not coexist: both benefit from intensity upgrades. I find it surprising, in 
view of the marvelous neutrino physics legacy of Los Alamos, that the community has not 
fought harder to keep the LAMPF/LANSCE neutrino program alive. 

Before closing, I would like to mention some other recent work relevant to LSND where 
neutron and neutrino physics are intertwined. In nuclear astrophysics it is known that 
approximately half of the elements heavier than Fe must be synthesized in the rapid-neutron- 
capture process, or r-process. The conditions for this process are very exotic, requiring 
explosive neutron fiuxes on the order of 10^^/cm^s lasting for ~ 1 s [10]. Under these 
conditions neutron capture becomes faster than /3 decay, so that the synthesis follows a 
path through very neutron rich nuclei defined by {n, 7) <-> (7, n) equilibrium. When the 
neutron flux diminishes, these isotopes then decay back to the valley of stability, forming 
the r-process nuclei we identify in nature. 

The question of the site of the r-process has been debated for almost 40 years. One of 
the most popular suggestions has been in the neutronized matter near the mass cut that is 
blown off the protoneutron star during a core-collapse supernova. This is the site relevant 



The protoneutron star cools by emitting neutrinos of all flavors. While the neutrinos 
remain in flavor equilibrium (due to neutral current interactions) throughtout most of their 
random walk out of the dense neutron star core, as they reach the neutrino sphere (the radius 
of last scattering) at ~ 10^^ g/cm^, some differences arise. The electron-flavor neutrinos are 
more strongly coupled to the matter because of the charged current interactions 



to LSND. 



1^^ + n 



p + e 



(7) 



n + e 



,+ 



(8) 
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As a result, the distributions of heavy- flavor neutrinos are characterized by temperatures 
T^^ ~ T^^ ~ 8 McV, while T^-^ - 5 MeV and T^^ ~ 4 MeV. (The latter difference arises 
because the matter is neutron rich.) 

Now if neutrinos have a cosmologically interesting mass (m^,^ or tjIi,^ ~ 2—50 eV) and even 
very modest mixing angles (sin^26' > 10~^), neutrinos can destroy the supernova r-process. 
Matter-enhanced oscillations between, for example, u^s and u^-s cause a temperature reversal 
between these species. As Qian, Fuller, and collaborators have discussed [17], the hot z/gS 
then have substantially increased charge current interactions, driving the matter above the 
neutron star proton rich. The r-process could no longer take place. 

Much of the oscillation parameter space suggested by LSND would then be eliminated by 
the requirement of a successful r-process. Of course, there appears to be an escape clause: 
perhaps the hot bubble above the protoneutron star is not the site of the r-process. As many 
other sites have been suggested, this is plausible and would invalidate any constraints on 
LSND. 

Recently Langanke, Qian, Vogel and I [18] explored another aspect of this r-process site, 
the fact that the neutrino flux would continue for several seconds after the r-process has been 
completed (which occurs when the matter temperature drops to ~ lO^K). These late neutri- 
nos would occasionally interact, by charged or neutral currents, with the r-process products. 
Typical reactions produce highly excited nuclei that decay by emitting several neutrons. 
Thus this neutrino postprocessing could alter the r-process abundance distribution. 

When we studied the postprocessing in detail, we discovered that there exist 8 rare nuclei 
(Te, W, and Re isotopes) that are extremely sensitivity to the neutrino postprocessing. These 
isotopes live in the abundance valleys immediately below the r-process abimdance peaks 
associated with the closed neutron shells at N = 50 and 82. Thus a few neutrino interactions 
with the isotopes in the abundance peaks can produce appreciable quantities of Te, W, and 
Re by spallation. When this possibility was explored in detailed, it turned out that the 
abundances of all eight isotopes could be entirely accounted for by neutrino spallation: the 
agreement with the known abundances is excellent (~ Icr). Furthermore the neutrino fiuence 
required to produce these nuclei by spallation is in good accord with supernova models of 
the r-process. Some results are illustrated in Fig. 1. 

This result strongly argues that these isotopes are neutrino postprocessing products, and 
therefore that the r-process does take place in the intense neutrino flux in the atmosphere 
above a protoneutron star. It follows that the r-process limit on neutrino oscillations [17] 
then must be taken as a very serious additional constraint on the same parameter space 
probed by LSND. 

5 Closing 

Let me close this session by saying "Thanks!" to 

* LAMPF, for 20 great years 

* Louis Rosen, for making it happen 

* our workshop hosts, Peter Barnes and Vernon Hughes, for allowing us to come together 
one last time 

* the organizing committee Cy Hoffman (chair), Ben Gibson, Bill Louis, Susan Ramsey, 
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Susan Seestrom, and Dave Vieira, for putting together a wonderful two days. 
This work was supported in part by the U.S. Department of Energy. 



References 

1. R. L. Huffman et al, Phys. Lett. 139B, 249 and references therein. 

2. I. S. Towner, A Shell Model Description of Light Nuclei (Clarendon Press, Oxford, 
1977) 383p. 

3. N. M. Hugenholtz, Physica 23, 533 (1957). 

4. T. H. Schucan and H. A. WeidenmuUer, Ann. Phys. 73, 108 (1972) and 76, 483 
(1973). 

5. W. C. Haxton and C. Johnson, Phys. Rev. Lett. 65, 1325 (1990). 

6. C. Bloch and J. Horowitz, Nucl. Phys. 8, 91 (1958). 

7. R. Haydock, J. Phys. A7, 2120 (1974); W. C. Haxton, in Proc. Int. Conf. on Nuclear 
Structure at High Angular Momentum, vol. 1, p. 294 (1993). 

8. K. F. Smith, et al, Phys. Lett. B234, 191 (1990). 

9. I. S. Alterev et al, Phys. Lett. B276, 242 (1992). 

10. J. P. Jacobs, W. M. Klipstein, S. K. Lamoreaux, B. R. Heckel, and E. N. Fortson, 
Phys. Rev. Lett. 71, 3782 (1993). 

11. G. Feinberg, Trans. N.Y. Acad. Sci. 38, 6 (1977); W. C. Haxton and E. M. Henley 
Phys. Rev. Lett. 51, 1937 (1983); V. V. Flambaum, I. B. Khriplovich, and O. P. Sushkov, 
Sov. Phys. JETP 60, 873 (1984); Phys. Lett. 162B, 213 (1985); Nucl. Phys. A449, 750 
(1986). 

12. L. 1. Schiff, Phys. Rev. 132, 2194 (1963). 

13. Ref. [10]; K. Lee, Ph.D thesis, Univ. of Washington (1994). 

14. N. Fortson, Phys. Rev. Lett. 70, 2383 (1993). 

15. C. Athanassopoulos et al, Phys. Rev. Lett. 77, 3082 (1996). 

16. See, for example, G. J. Mathews and J. J. Cowan, Nature (London) 345, 491 (1990), 
and references therein. 

17. Y.-Z. Qian et al, Phys. Rev. Lett. 71, 1965 (1993); G. M. Fuller, Phys. Rep. 227, 
149 (1993). 

18. W. C. Haxton, K. Langanke, Y. -Z. Qian, and P. Vogel, Phys. Rev. Lett. 78, 2694 
(1997) and Phys. Rev. C. 55, 1532 (1997). 



8 



Figure Captions 

Figure 1: The solid line represents the r-process distribution that would be produced if, 
after freezeout, there is no further processing by neutrinos. The dotted line 
represents the changes due to the postprocessing by neutrinos (corresponding to 
a total neutrino fluence of 0.015 • 10^^ erg/km^). The fluence was chosen so that 
isotopes in the window A = 183-187 could be produced entirely by neutrinos. 
The quality of the resulting fit to abundances in this window strongly suggests 
that the r-process occurs in an intense neutrino flux. See Ref. [18] for more 
details. 
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